Sources of biases in the in vitro testing of nanomaterials: the role of the biomolecular corona.

The biological fate of nanomaterials (NMs) is driven by specific interactions through which biomolecules, naturally adhering onto their surface, engage with cell membrane receptors and intracellular organelles. The molecular composition of this layer, called the biomolecular corona (BMC), depends on both the physical-chemical features of the NMs and the biological media in which the NMs are dispersed and cells grow. In this work, we demonstrate that the widespread use of 10% fetal bovine serum in an in vitro assay cannot recapitulate the complexity of in vivo systemic administration, with NMs being transported by the blood. For this purpose, we undertook a comparative journey involving proteomics, lipidomics, high throughput multiparametric in vitro screening, and single molecular feature analysis to investigate the molecular details behind this in vivo/in vitro bias. Our work indirectly highlights the need to introduce novel, more physiological-like media closer in composition to human plasma to produce realistic in vitro screening data for NMs. We also aim to set the basis to reduce this in vitro-in vivo mismatch, which currently limits the formulation of NMs for clinical settings.

Molecular Dynamics and Machine Learning Give Insights on the Flexibility-Activity Relationships in Tyrosine Kinome.

Tyrosine kinases are a subfamily of kinases with critical roles in cellular machinery. Dysregulation of their active or inactive forms is associated with diseases like cancer. This study aimed to holistically understand their flexibility-activity relationships, focusing on pockets and fluctuations. We studied 43 different tyrosine kinases by collecting 120 µs of molecular dynamics simulations, pocket and residue fluctuation analysis, and a complementary machine learning approach. We found that the inactive forms often have increased flexibility, particularly at the DFG motif level. Noteworthy, thanks to these long simulations combined with a decision tree, we identified a semiquantitative fluctuation threshold of the DGF+3 residue over which the kinase has a higher probability to be in the inactive form.

On the allosteric puzzle and pocket crosstalk through computational means.

Allostery is a constitutive, albeit often elusive, feature of biomolecular systems, which heavily determines their functioning. Its mechanical, entropic, long-range, ligand, and environment-dependent nature creates far from trivial interplays between residues and, in general, the secondary structure of proteins. This intricate scenario is mirrored in computational terms as different notions of "correlation" among residues and pockets can lead to different conclusions and outcomes. In this article, we put on a common ground and challenge three computational approaches for the correlation estimation task and apply them to three diverse targets of pharmaceutical interest: the androgen A2A receptor, the androgen receptor, and the EGFR kinase domain. Results show that partial results consensus can be attained, yet different notions lead to pointing the attention to different pockets and communications.

Probabilistic Pocket Druggability Prediction via One-Class Learning.

The choice of target pocket is a key step in a drug discovery campaign. This step can be supported by in silico druggability prediction. In the literature, druggability prediction is often approached as a two-class classification task that distinguishes between druggable and non-druggable (or less druggable) pockets (or voxels). Apart from obvious cases, however, the non-druggable class is conceptually ambiguous. This is because any pocket (or target) is only non-druggable until a drug is found for it. It is therefore more appropriate to adopt a one-class approach, which uses only unambiguous information, namely, druggable pockets. Here, we propose using the import vector domain description (IVDD) algorithm to support this task. IVDD is a one-class probabilistic kernel machine that we previously introduced. To feed the algorithm, we use customized DrugPred descriptors computed via NanoShaper. Our results demonstrate the feasibility and effectiveness of the approach. In particular, we can remove or mitigate biases chiefly due to the labeling.

Probing Interplays between Human XBP1u Translational Arrest Peptide and 80S Ribosome.

The ribosome stalling mechanism is a crucial biological process, yet its atomistic underpinning is still elusive. In this framework, the human XBP1u translational arrest peptide (AP) plays a central role in regulating the unfolded protein response (UPR) in eukaryotic cells. Here, we report multimicrosecond all-atom molecular dynamics simulations designed to probe the interactions between the XBP1u AP and the mammalian ribosome exit tunnel, both for the wild type AP and for four mutant variants of different arrest potencies. Enhanced sampling simulations allow investigating the AP release process of the different variants, shedding light on this complex mechanism. The present outcomes are in qualitative/quantitative agreement with available experimental data. In conclusion, we provide an unprecedented atomistic picture of this biological process and clear-cut insights into the key AP-ribosome interactions.

Opportunities and Challenges for Machine Learning in Rare Diseases.

Rare diseases (RDs) are complicated health conditions that are difficult to be managed at several levels. The scarcity of available data chiefly determines an intricate scenario even for experts and specialized clinicians, which in turn leads to the so called "diagnostic odyssey" for the patient. This situation calls for innovative solutions to support the decision process via quantitative and automated tools. Machine learning brings to the stage a wealth of powerful inference methods; however, matching the health conditions with advanced statistical techniques raises methodological, technological, and even ethical issues. In this contribution, we critically point to the specificities of the dialog of rare diseases with machine learning techniques concentrating on the key steps and challenges that may hamper or create actionable knowledge and value for the patient together with some on-field methodological suggestions and considerations.

Machine Learning and Enhanced Sampling Simulations for Computing the Potential of Mean Force and Standard Binding Free Energy.

Computational capabilities are rapidly increasing, primarily because of the availability of GPU-based architectures. This creates unprecedented simulative possibilities for the systematic and robust computation of thermodynamic observables, including the free energy of a drug binding to a target. In contrast to calculations of relative binding free energy, which are nowadays widely exploited for drug discovery, we here push the boundary of computing the binding free energy and the potential of mean force. We introduce a novel protocol that leverages enhanced sampling, machine learning, and ad hoc algorithms to limit human intervention, computing time, and free parameters in free energy calculations. We first validate the method on a host-guest system, and then we apply the protocol to glycogen synthase kinase 3 beta, a protein kinase of pharmacological interest. Overall, we obtain a good correlation with experimental values in relative and absolute terms. While we focus on protein-ligand binding, the strategy is of broad applicability to any complex event that can be described with a path collective variable. We systematically discuss key details that influence the final result. The parameters and simulation settings are available at PLUMED-NEST to allow full reproducibility.

Molecular Recognition by Gold Nanoparticle-Based Receptors as Defined through Surface Morphology and Pockets Fingerprint.

Ligand shell-protected gold nanoparticles can form nanoreceptors that recognize and bind to specific molecules in solution, with numerous potential innovative applications in science and industry. At this stage, the challenge is to rationally design such nanoreceptors to optimize their performance and boost their further development. Toward this aim, we have developed a new computational tool, Nanotron. This allows the analysis of molecular dynamics simulations of ligand shell-protected nanoparticles to define their exact surface morphology and pocket fingerprints of binding cavities in the coating monolayer. Importantly, from dissecting the well-characterized pairing formed by the guest salicylate molecule and specific host nanoreceptors, our work reveals that guest binding at such nanoreceptors occurs via preformed deep pockets in the host. Upon the interaction with the guest, such pockets undergo an induced-fit-like structural optimization for best host-guest fitting. Our findings and methodological advancement will accelerate the rational design of new-generation nanoreceptors.

Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation.

Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.

Spathial: an R package for the evolutionary analysis of biological data.

SUMMARY: A primary problem in high-throughput genomics experiments is finding the most important genes involved in biological processes (e.g. tumor progression). In this applications note, we introduce spathial, an R package for navigating high-dimensional data spaces. spathial implements the Principal Path algorithm, which is a topological method for locally navigating on the data manifold. The package, together with the core algorithm, provides several high-level functions for interpreting the results. One of the analyses we propose is the extraction of the genes that are mainly involved in the progress from one state to another. We show a possible application in the context of tumor progression using RNA-Seq and single-cell datasets, and we compare our results with two commonly used algorithms, edgeR and monocle3, respectively. AVAILABILITY AND IMPLEMENTATION: The R package spathial is available on the Comprehensive R Archive Network (https://cran.r-project.org/web/packages/spathial/index.html) and on GitHub (https://github.com/erikagardini/spathial). It is distributed under the GNU General Public License (version 3). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multi-target dopamine D3 receptor modulators: Actionable knowledge for drug design from molecular dynamics and machine learning.

Local changes in the structure of G-protein coupled receptors (GPCR) binders largely affect their pharmacological profile. While the sought efficacy can be empirically obtained by introducing local modifications, the underlining structural explanation can remain elusive. Here, molecular dynamics (MD) simulations of the eticlopride-bound inactive state of the Dopamine D3 Receptor (D3DR) have been clustered using a machine learning-based approach in the attempt to rationalize the efficacy change in four congeneric modulators. Accumulating extended MD trajectories of receptor-ligand complexes, we observed how the increase in ligand flexibility progressively destabilized the crystal structure of the inactivated receptor. To prospectively validate this model, a partial agonist was rationally designed based on structural insights and computational modeling, and eventually synthesized and tested. Results turned out to be in line with the predictions. This case study suggests that the investigation of ligand flexibility in the framework of extended MD simulations can assist and inform drug design strategies, highlighting its potential role as a powerful in silico counterpart to functional assays.

Investigating Drug-Target Residence Time in Kinases through Enhanced Sampling Simulations.

It is widely accepted that drug-target association and dissociation rates directly affect drug efficacy and safety. To rationally optimize drug binding kinetics, one must know the atomic arrangement of the protein-ligand complex during the binding/unbinding process in order to detect stable and metastable states. Whereas experimental approaches can determine kinetic constants with fairly good accuracy, computational approaches based on molecular dynamics (MD) simulations can deliver the atomistic details of the unbinding process. Furthermore, they can also be utilized prospectively to predict residence time (i.e., the inverse of unbinding kinetics constant, koff) with an acceptable level of accuracy. Here, we report a novel method based on adiabatic bias MD with an electrostatics-like collective variable (dubbed elABMD) for sampling protein-ligand dissociation events in two kinases. elABMD correctly ranked a ligand series on glucokinase, in agreement with experimental data and previous calculations. Subsequently, we applied the new method prospectively to a congeneric series of GSK-3ß inhibitors. For this series, new crystal structures were generated and the residence time was experimentally measured with surface plasmon resonance (SPR). There was good agreement between computational predictions and experimental measures, suggesting that elABMD is an innovative and efficient tool for calculating residence times.

Finding Principal Paths in Data Space.

In this paper, we introduce the concept of principal paths in data space; we show that this is a well-characterized problem from the point of view of cognition, and that it can lead to salient insights in the analyzed data enabling topological/holistic descriptions. These paths, interestingly, can be interpreted as local principal curves, and in this paper, we suggest that they are analogous to what, in the statistical mechanics realm, are called minimum free-energy paths. Here, we move that concept from physics to data space and compute them in both the original and the kernel space. The algorithm is a regularized version of the well-known k -means clustering algorithm. The regularization parameter is derived via an in-sample model selection process based on the Bayesian evidence maximization. Interestingly, we show that this choice for the regularization parameter consistently leads to the same manifold even when changing the number of clusters. We apply the method to common data sets, dynamical systems, and, in particular, to molecular dynamics trajectories showing the generality, the usefulness of the approach and its superiority with respect to other related approaches.

NanoShaper-VMD interface: computing and visualizing surfaces, pockets and channels in molecular systems.

SUMMARY: NanoShaper is a program specifically aiming the construction and analysis of the molecular surface of nanoscopic systems. It uses ray-casting for parallelism and it performs analytical computations whenever possible to maximize robustness and accuracy of the approach. Among the other features, NanoShaper provides volume, surface area, including that of internal cavities, for any considered molecular system. It identifies pockets via a very intuitive definition based on the concept of probe radius, intrinsic to the definition of the solvent excluded surface. We show here that, with a suitable choice of the parameters, the same approach can also permit the visualisation of molecular channels. NanoShaper has now been interfaced with the widely used molecular visualization software VMD, further enriching its already well furnished toolset. AVAILABILITY AND IMPLEMENTATION: VMD is available at http://www.ks.uiuc.edu/Research/vmd/. NanoShaper, its documentation, tutorials and supporting programs are available at http://concept.iit.it/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BiKi Life Sciences: A New Suite for Molecular Dynamics and Related Methods in Drug Discovery.

In this paper, we introduce the BiKi Life Sciences suite. This software makes it easy for computational medicinal chemists to run ad hoc molecular dynamics protocols in a novel and task-oriented environment; as a notebook, BiKi (acronym of Binding Kinetics) keeps memory of any activity together with dependencies among them. It offers unique accelerated protein-ligand binding/unbinding methods and other useful tools to gain actionable knowledge from molecular dynamics simulations and to simplify the drug discovery process.

Fast Dynamic Docking Guided by Adaptive Electrostatic Bias: The MD-Binding Approach.

Engineering chemical entities to modify how pharmaceutical targets function, as it is done in drug design, requires a good understanding of molecular recognition and binding. In this context, the limitations of statically describing bimolecular recognition, as done in docking/scoring, call for insightful and efficient dynamical investigations. On the experimental side, the characterization of dynamical binding processes is still in its infancy. Thus, computer simulations, particularly molecular dynamics (MD), are compelled to play a prominent role, allowing a deeper comprehension of the binding process and its causes and thus a more informed compound selection, making more significant the computational contribution to drug discovery (Carlson, H. A. Curr. Opin. Chem. Biol. 2002, 6, 447-452). Unfortunately, MD-based approaches cannot yet describe complex events without incurring prohibitive time and computational costs. Here, we present a new method for fully and dynamically simulating drug-target-complex formations, tested against a real world and pharmaceutically relevant benchmark set. The method, based on an adaptive, electrostatics-inspired bias, envisions a campaign of trivially parallel short MD simulations and a strategy to identify a near native binding pose from the sampled configurations. At an affordable computational cost, this method provided predictions of good accuracy also when the starting protein conformation was different from that of the crystal complex, a known hurdle for traditional molecular docking (Lexa, K. W.; Carlson, H. A. Q. Rev. Biophys. 2012, 45, 301-343). Moreover, along the observed binding routes, it identified some key features also found by much more computationally expensive plain-MD simulations. Overall, this methodology represents significant progress in the description of binding phenomena.

Allosteric Communication Networks in Proteins Revealed through Pocket Crosstalk Analysis.

The detection and characterization of binding pockets and allosteric communication in proteins is crucial for studying biological regulation and performing drug design. Nowadays, ever-longer molecular dynamics (MD) simulations are routinely used to investigate the spatiotemporal evolution of proteins. Yet, there is no computational tool that can automatically detect all the pockets and potential allosteric communication networks along these extended MD simulations. Here, we use a novel and fully automated algorithm that examines pocket formation, dynamics, and allosteric communication embedded in microsecond-long MD simulations of three pharmaceutically relevant proteins, namely, PNP, A2A, and Abl kinase. This dynamic analysis uses pocket crosstalk, defined as the temporal exchange of atoms between adjacent pockets, along the MD trajectories as a fingerprint of hidden allosteric communication networks. Importantly, this study indicates that dynamic pocket crosstalk analysis provides new mechanistic understandings on allosteric communication networks, enriching the available experimental data. Thus, our results suggest the prospective use of this unprecedented dynamic analysis to characterize transient binding pockets for structure-based drug design.

Import Vector Domain Description: A Kernel Logistic One-Class Learning Algorithm.

Recognizing the samples belonging to one class in a heterogeneous data set is a very interesting but tough machine learning task. Some samples of the data set can be actual outliers or members of other classes for which training examples are lacking. In contrast to other kernel approaches present in the literature, in this work, the problem is faced defining a one-class kernel machine that delivers the probability for a sample to belong to the support of the distribution and that can be efficiently trained by a hybrid sequential minimal optimization-expectation maximization algorithm. Due to the analogy to the import vector machine and to the one-class approach, we named the method import vector domain description (IVDD). IVDD was tested on a toy 2-D data set in order to characterize its behavior on a set of widely used benchmarking UCI data sets and, lastly, challenged against a real world outlier detection data set. All the results were compared against state-of-the-art closely related methods such as one-class-SVM and Support Vector Domain Description, proving that the algorithm is equally accurate with the additional advantage of delivering the probability estimate for each sample. Finally, a few variants aimed at providing memory savings and/or computational speed-up in the light of big data analysis are briefly sketched.